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A hybrid orbital model for the topological insulator Bi2Se3 is proposed and compared to that of 
graphene. The existence of a 7r-bond trimer on the Sel layer at the van der Waals gap and on the 
surface is thought to be responsible for the unusual surface conduction mechanism found in Bi2Se3 
as a topological insulator. The three 7r-bonds are locked as permanent electric dipoles between the 
Sel layers in the van der Waals gap as an attractive force for the bulk as a gapped semiconductor. 
In contrast, the 7r-bond trimers on the surface are unlocked and exhibit two degenerate quantum 
states, creating a local charge transfer mechanism in the cross-bridge model that is responsible for 
the surface conduction. The role of 7r-bonds on the Bi2Se3 surface is compared with that in graphene 
for a similar 2D band structure containing Dirac cones. 

PACS numbers: 73.20.-r; 73.25.+i; 31.10.+z; 33.15.Fm 



I. INTRODUCTION 

Topological insulator is a unique state of matter for a 
band insulator with surface-only conduction and a spin- 
polarized surface state. Bi2Se3 is a the second-generation 
topological insulator with a nearly idealized single Dirac 
cone at the T-pointi The crystal structure of Bi2Se3 
(Fig. [1]) can be viewed from the Bi-Se quintuple layers 
of hep close packing in space group of R3m with van der 
Waals gaps between the Sel-layers^ The observation of 
gapless Dirac cone in the surface band structure of topo- 
logical insulators and in the perfect 2D electron gas of 
graphene implies that a similar conduction mechanism is 
expected^ However, the understanding of this new class 
of material has been based on analogous descriptions of 
topological invariants, Berry's phase, and the quantum 
spin Hall effect. We believe a real space view from the 
chemical bond perspective should enhance and compli- 
ment to the material design significantly. In particular, 
we find that there is a great deal of similarity between 
the graphene and Bi2Se3 surface in terms of the Dirac 
cone band picture and the existence of 7r-bonds. 



II. CURRENT VIEW OF CRYSTAL AND 
ELECTRONIC STRUCTURES FOR Bi 2 Se 3 AND 
GRAPHENE 



Graphite has three neighboring carbon atoms 
([He]2s 2 2p 2 ) in each layer, with graphene defined as a 
monolayer of graphite, which leads to the hybridization 
of triangular-shaped sp 2 orbitals with three coordination 
sites with the three neighboring carbons in the a-bonds 
and leaves one unpaired electron in the p-orbital to form 
7r-bond with one of the three neighboring carbon atoms 
statistically^ Whereas 7r-bonds can be formed with any 
one of the three equivalent neighboring carbon atoms, 
instead of choosing a particular one to bond with, the 



bond order is categorized as l| among all carbon atoms 
on a time-average base£ The addition and the dynamic 
nature of 7r-bonds cause graphene to show a stronger 
bond strength that is manifested by a higher melting 
point compared with that of diamond; diamond only has 
one a-bond among the carbon atoms in the tetrahedral- 
shaped sp 3 orbitals and can be categorized as having 
a bond order of 1. The band picture of graphene has 
often been described by the uniquely connected Dirac 
cones between the valence and conduction bands, in 
which linear dispersion of E(k) suggests the conduction 
mechanism is describable by the massless Dirac fermions 
with a relativistic mass of ~ 0.007 m e and a Fermi 
velocity of ~ c/300. 6r8 The conduction of graphite 
has often been interpreted through the electron carrier 
transport in the anti-bonding 7r*-band as the conduction 
band. 15 Perhaps a more appropriate description of the 
conduction in graphene would be a chiral tunneling 
mechanism, which has also been suggested to be the first 
example of the "Klein paradox" in condensed matter^ 

The interpretation of how can a topological insulator 
surface conduct electricity relies heavily on the argument 
of the band inversion phenomenon proposed in heavy el- 
ements with strong spin-orbit coupling^ band inversion 
for the partially filled p orbitals in Bi and Se may exist 
under the crystal field field splitting and strong spin-orbit 
coupling. Considering the medium electronegativity dif- 
ference of ~0.5 between Bi and Se, orbital hybridiza- 
tion to form polar covalent bonds among Bi-Se is possi- 
ble. The hybrid orbital model for semiconductor silicone 
has been fully accepted to have four sp 3 hybridized a- 
bonds in tetrahedral coordination among the Si atoms of 
atomic electronic configuration of [Ne]3s 2 3p 2 , where sili- 
cone crystal has a band gap of ~1 eV, just like carbon in 
the diamond structure. Since Bi2Se3 has a band gap size 
of ~0.3 eV in the semiconductor range and both Bi/Se 
elements are in close packing with octahedral coordina- 
tion, we believe that there is no reason not to describe 
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FIG. 1. (Color online) The crystal structure of Bi2Se3 is 
shown with two quintuple layers along the [110] direction for 
simplicity. The insets show the actual octahedral environ- 
ment for Sel at the van der Waals gap, Se2 in the middle, 
and Bi from left to right, respectively. 



Bi2Se3 with fully hybridized s-p-d orbitals, especially if 
the newly proposed hybridization is able to explain the 
defining physical properties of topological insulators. In 
particular, we find the existence of surface 7r-bond is un- 
avoidable under the newly proposed sp 3 d 2 hybridization, 
and 7r-bond is the key component responsible for the elec- 
tronic structure and physical property interpretation for 
Bi2Se3 as a topological insulator, which shows great sim- 
ilarity to that of graphene. 



III. CRYSTAL STRUCTURE AND CHEMICAL 
BONDS OF Bi 2 Se 3 

The crystal structure of Bi2Se3 can be described best 
using a Bi-Se quintuple as a layer unit, and these quin- 
tuple layers are close-packed with symmetry describable 
with the R3m space group, i.e., three sliding planes per 
unit c-axis, as shown in Fig. [TJ Herein, we propose a 
real space view of Bi2Se3 through chemical bonds con- 
structed in the hybrid orbital model and use it to explain 
the defining physical properties observed in Bi2Se3 as a 
topological insulator. We will tackle the special physical 
properties found in Bi2Se3 through real space chemical 
bonding instead of using the band structure approach. 
Parallel interpretations in both real and reciprocal spaces 
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FIG. 2. (Color online) The proposed orbital hybridization for 
one quintuple layer of Bi2Se3- The green dashed lines rep- 
resent the six hybridized channels of sp 3 d 2 in the octahedral 
shape, in which cr-bonds are formed between Bi-Se, but the 
three unpaired 7r-electrons for Sel at the van der Waals gap 
or surface may form a 7r-bond trimer with three out of six Sel 
atoms in the same layer. 



should greatly enhance the understanding of topological 
insulators and lead to the identification of more materials 
in this class. 

Considering the number of outer-shell electrons in 
Bi([Xe]4f 14 5d 10 6s 2 6p 3 ) and Se([Ar]3d 10 4s 2 4p 4 ) and the 
close-packed BiSe6-SeBi6 octahedra, the chemical bond- 
ing of Bi2Se3 has been proposed based on different choices 
of orbital hybridization before, such as p 3 d 3 , sp 3 d 2 , or a 
combination of the two i 1Q i n Early proposed hybridiza- 
tion models described the unit quintuple layer only with- 
out considering the relatively strong van der Waals force. 
In addition, crucial information about the Sel layer sur- 
face is missing. We anticipate that a correct hybrid or- 
bital model for Bi2Se3 should be able to explain the rel- 
atively strong van der Waals force, the bulk as a gapped 
semiconductor, and the unique surface conduction as a 
topological insulator. 

In each quintuple layer unit of Bi2Se3, there are two 
different symmetry sites for the Se atoms, one is in the 
center (Se2), and the other two are at the van der Waals 
gap (Sel), as shown in Fig. [TJ Although there are six 
atoms surround the central one in octahedral coordina- 
tion for both Bi and Se, the available s-p-d orbitals for 
both the Bi and Se atoms should allow sp 3 d 2 hybridiza- 
tion, as shown in Fig. [2j The six outer-shell electrons in 
Se can easily satisfy the requirement of the six sp 3 d 2 hy- 
bridized orbitals, but we note that the outer-shell of Bi 
has only five electrons. One of the 5d 10 electrons should 
be allowed to move up from the filled bd orbitals to par- 
ticipate in the sp 3 d 2 hybridization naturally in exchange 
for the bonding energy gain, especially when the energy 
difference between the n = 5 and n = 6 levels decreases, 
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FIG. 3. (Color online) The conjugated systems formed with tt- 
bonds for (a) graphene, (b) Sel layer on the surface of Bi2Se3, 
and (c) benzene (CeHe). Graphene has one 7r-bond shared 
by three carbon atoms in cr-bonds (sp 2 ) as illustrated on the 
right. The Sel layer on the surface of Bi2Se3 has three tt- 
bonds, forming a trimer with three of the six neighboring 
Se atoms. There are two degenerated states of 7r-bond trimer 
rotated 60° relative to each other. Benzene (CeHe) also shows 
two configurations with three 7r-bonds, but this is a closed 
conjugated system. 



as expected for heavy elements with a large principal 
quantum number n. It is clear that both the Bi and Se2 
can satisfy the six cr-bonds required for octahedral shape 
coordination and create one well-defined quintuple layer. 
However, we note that the Sel at the van der Waals gap 
will have three electrons remain unpaired. The three un- 
paired electrons on Sel surface should separate as far as 
possible to 120° apart in a plane to form a trimer to 
reduce the Coulumb repulsion between them, and three 
7r-bonds can be formed with three out of the six neigh- 
boring Sel atoms in the same layer. We call the three 
7r-bonds in one unit as a 7r-bond trimer, and a 7r-electron 
is used to denote the unpaired electron after the 7r-bond 
is broken. 



IV. tt-BOND CONJUGATED SYSTEMS 

For graphene, it is well established in chemistry that 
statistically, only one unpaired 7r-electron in the half- 
filled p^-orbital would form a 7r-bond with one of three 
neighboring carbon atoms and overlap with one of the 
three cr-bonds from the hybridized sp 2 orbit als, as shown 
in Fig. [3{a)ji£ The derealization of the 7r-electron af- 
ter 7r-bond breaking would form a 7r-bond with another 
neighboring carbon instantaneously, and the three tt- 
bond bonding/nonbonding choices form one conjugated 
system. Although it seems possible for all 7r-bonds to 
orderly lock when considering enthalpy reduction only, 
this would be less preferable considering the lowering 
of the total Gibbs free energy for the system that con- 
tains competing enthalpy and entropy terms, i.e., the tt- 
electron configurational entropy gain must also be con- 
sidered for the system containing easily perturbed, ex- 
tremely small energy gain/loss through the 7r-bond bond- 
ing/nonbonding energy exchange. For Bi2Se3, the Sel 
surface can be exposed after cleaving the crystal through 
the van der Waals gap. After the Sel layer is exposed, 
the three unpaired 7r-electrons will form a 7r-bond trimer 
resonating between the two equivalent states as one con- 
jugated system, as illustrated in Fig. [3fb). 

The entropy factor must be taken into account to de- 
termine the comparable impact for a system searching for 
the true ground state. In the language of chemical bond- 
ing, electrons localized in the valence band are locked in 
cr-bonds, and the semi- localized electrons on the Bi2Se3 
surface correspond to the much weaker 7r-bond. The tt- 
bond is constructed through side-to-side minimum or- 
bital overlapping (see Fig. [3j) and is so weak that sponta- 
neous bond breaking is allowed^ Another good example 
is benzene (CeHe), which has three 7r-bonds shared by 
six carbon atoms. The two equivalent choices of tt— bond 
configurations and the small 7r-bond bonding energy do 
not preferentially cause the system to settle into either 
one of the two arrangements, as shown in Fig. [3jc)^ 
Instead, a dynamic conjugated system of raised configu- 
rational entropy is the best choice. 

We might describe such a dynamic 7r-bond energy ex- 
change in time and space using a Feynman diagram-like 
picture, as illustrated in Fig.HJ Similar entropy-enthalpy 
competition has been applied to the uncertain 7r-bond po- 
sition among the three choices of the directional cr-bonds, 
as shown in Fig. 02a) for graphene, i.e., although it takes 
energy to break the existing 7r-bond at a specific time 
and position, the new 7r-bond formation at a different 
position with a configurational entropy gain would still 
be able to lower the total Gibbs free energy for the re- 
maining terms within AG = AH — TAS. In the language 
of thermodynamics, which addresses the initial and final 
states only, AH=0 always, although the 7r-bond bond- 
ing energy is nonzero and exchanged constantly. In na- 
ture, the dynamic conjugated system is preferred over 
the static choices (whether ordered or not) for systems 
with energy gaps as small as the weak tt— bond bonding 
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FIG. 4. (Color online) An illustration of the time and space 
relationship for 7r-bond, 7r-electron, and 7r-bond bonding en- 
ergy exchange as a quasiparticle with a Feynman diagram-like 
picture. 



Fig. [TJ the bond length of Sel-Bi (2.85A) in the virtual 
octahedron of (Sel)Bis(Sel)3 is significantly different to 
that of Se2-Bi (3.08A) in the (Se2)Bi 6 octahedron^ In 
addition, the Sel atoms are arranged in a 2D triangu- 
lar lattice, i.e., there are six Sel atoms surrounding each 
Sel, as shown on the opposite sides of the van der Waals 
gap shown in Fig. [5] For each Sel atom exposed to the 
van der Waals gap, there are three unpaired electrons 
waiting to bond. These electrons are called 7r-electrons 
before a 7r-bond is formed with the neighboring Sel. The 
three half-filled sp 3 <i 2 -orbitals can form three 7r-bonds 
with three of the six neighboring Sel atoms on the sur- 
face, with a weaker bonding energy due to the minimum 
orbital overlapping compared with the a-bond. Consid- 
ering the Coulomb repulsion among the three unpaired 
electrons, a symmetric triangular distribution should be 
preferred with two available degenerate states, which is 
defined as a 7r-bond trimer for simplicity. When two Sel 
layers approach each other in the van der Waals gap, a 
total of six electric dipoles can form in such a way that 
two 7r-bond trimers are locked in a rotation that is 60° 
relative to the other, which could be viewed as the source 
of an attractive force in the van der Waals gap. 




FIG. 5. (Color online) The electric dipoles formed with tt- 
bonds in the van der Waals gap between (a) two Sel layers in 
the van der Waals gap of Bi2Se3, and (b) between the bilayer 
graphene. All six electric dipoles are locked definitely for the 
former, whereas one of the three possible dipoles must be 
frustrated in bilayer graphene to weaken its attractive force. 
Note that these two layers are actually shifted for hexagonal 
close packing in reality. 



energy. 



V. VAN DER WAALS GAP IN Bi 2 Se 3 AND 
BILAYER GRAPHENE 

Although the quintuple layers are stacked with sliding 
2D triangular lattices to form close packing, as shown in 



It is interesting to compare the strength of the van 
der Waals force between the bilayer graphene and Bi2Se3 
through 7r-bond and electric dipole formation. Consid- 
ering that the surface of Bi2Se3 is composed of com- 
pletely free 7r-bond trimers, each 7r-bond trimer has two 
possible choices in the next move through clockwise or 
counterclockwise motion. When a second identical Se 
surface approaches just as in van der Waals gap forma- 
tion, the second Se surface has no choice but to form 
six permanent electric dipoles to avoid Coulomb repul- 
sion and locks to the first Se surface accordingly, as 
shown in Fig. E^a). In contrast, the 7r-bond in each 
graphene layer has three choices at any instance. The 
approaching second graphene layer would always have 
one of the three possible dipoles frustrated relative to 
that of the first layer, which leads to a weaker attractive 
force. The weaker van der Waals force is the key prop- 
erty that makes the mechanical extraction of a monolayer 
of graphite (graphene) possible^ whereas we find that 
the same mechanical cleavage method failed to prepare 
a mono quintuple layer from single crystal Bi2Se3. Con- 
sidering the quintuple layer thickness of ~7 A and the 
van der Waals gap size of ~2.5 A, it is not surprising to 
find that for Bi2Se3, a gapless Dirac cone can be observed 
and defined as a topological insulator for only up to five 
layers of quintuples, 14 which must be closely related to 
the inter 7r-bond trimer distance across the actually ex- 
posed Sel surface, i.e., whether 7r-bond trimer is locked 
is closely related to the gap opening at the Dirac cone 
for surface conduction. We will describe this concept in 
greater detail in the following section. 
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FIG. 6. (Color online) A cross-bridge model can be used to 
explain the 2D conduction mechanism for a snapshot of (a) 
graphene, in which an (lcr+|7r)-bond is found among all car- 
bons on time- average basis, and (b) a Sel surface of Bi2Se3, 
in which the solid and empty lobes represent the two choices 
of 7r-bond trimers that lead to |7r-bonds on a time- average 
basis among all Sel atoms. An infinite percolation route can 
be constructed for three-way connection using one constantly 
rotating 7r-bond bridge for graphene, and for a six-way con- 
nection with one rotating 7r-bond trimer for Sel layers on the 
surface of Bi 2 Se3. 



VI. tt-BOND CONDUCTION MECHANISM 

We find it is necessary to re-examine the conduction 
mechanism of graphene and Bi2Se3 surfaces using the 
unique quantum nature of 7r-bonds in real space, instead 
of following the conventional view of the 7r* conduction 
band picture that treats unpaired p z electrons as itin- 
erant carriers. 15 Unlike the conventional view that dis- 
tinguishes electrons in either the localized state in the 
valence band or in the itinerant state in the conduction 
band, the 7r-electron has a unique dual role of coexist- 
ing localized (spin) and itinerant (charge) characteris- 
tics in the time window, as expected for the chiral elec- 
trons described by the quantum spin Hall effect, and the 
surface bandgap is so small that quantum tunneling be- 
comes possible, 7 which is also describable by the impact 
of entropy on the system. The constant 7r-bond bond- 
ing/nonbonding action among the three choices in space 
for each atom site locally can be illustrated with a Feyn- 
man diagram-like picture in space and time as shown in 
Fig. SI where 7r-bond bonding energy can be exchanged 
to form a quasiparticle-like conjugated system. 

The Dirac cone with a linear dispersion of E(k) also im- 
plies that the surface 2D conduction mechanism is better 
described by Dirac fermions following relativistic dynam- 
ics with zero effective mass* 7 - Such massless transport im- 
plies that itinerant electrons are not the actual carriers 



for electricity transport; instead, the spontaneous and 
constant energy exchange should be responsible for the 
conduction of electricity. To construct a continuous route 
for the electricity flow based on the local 7r-bond constant 
bonding/nonbonding action, as described in Fig. [4j we 
propose a cross-bridge model to explain the unique con- 
duction mechanism for both graphene and Bi2Se3 sur- 
faces in real space. 

In the case of graphene, only one 7r-bond can form 
with one of the three neighboring carbon atoms at any 
instance, but a continuous percolation route can always 
be formed at any instance, as shown in Fig. [6^a). In 
the case of the Sel layer on a Bi2Se3 surface, there is 
one 7r-bond trimer (three 7r-bonds) shared with the six 
neighboring Se atoms, and a continuous percolation route 
can also always form, as shown in Fig. [6f b). Graphene 
conduction can be described with one constantly shifting 
7r-bond bridge among the three a-bonds, which seems to 
have an (lcr+|7r)-bond between any two carbons on a 
time- average basis, but the action of 7r-bond shifting is 
actually discrete in time and localized on the same atom. 
Similarly, the bond order among Sel atoms on the sur- 
face of Bi2Se3 can also be viewed as a 2~^" - 

bond on a 

time- average basis, although the 7r-bond trimer is con- 
stantly shifting between the two choices of clockwise and 
counterclockwise rotation discretely in time, centered on 
each Sel atom. 

It is also interesting to note that such a cross-bridge 
model has been applied in the field of biomechanics to 
explain the continuous sliding function of muscles. The 
cross-bridge model explains how a myosin filament with 
only three myosin heads (separated 120°) can contin- 
uously attach the actin filaments arranged in a close- 
packed hexagonal 2D lattice in a cross-sectional view^ 
It must be noted that the existence of 7r-bonds does 
not guarantee electric conduction. For example, benzene 
(CqHg) has three 7r-bonds within one loop, as shown in 
Fig.[3fc). The two choices of 7r-bond configurations form 
one dynamic conjugated system of two equivalent states, 
but the open and close actions of 7r-bonds do not cre- 
ate an open and continuous percolation route for charge 
transfer^ In contrast, conduction becomes possible in 
a conducting polymer that contains both a 7r-bond con- 
jugated system and an open and continuous percolation 
routed 7 - 

How do graphene and the Bi2Se3 surface conduct elec- 
tricity in the absence of the conventionally defined itin- 
erant electron carriers? First, there is a percolation net- 
work of nonbonding 7r-electrons that dynamically exists 
all the time according to the cross-bridge model discussed 
above. In other words, the 7r-bond trimers must be un- 
locked from the electric dipole formed in the van der 
Waals gap first. Second, although the 7r-electron remains 
in the local vicinity of each atom, the energy is allowed to 
be transferred between these atoms through a relay-type 
energy exchange, i.e, although it takes 7r-bond bonding 
energy to excite the electron to the nearly-free 7r-electron 
state before it bonds again statistically, the same amount 
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of 7r-bonding energy is exchanged between atoms like a 
quasiparticle, as illustrated in Fig. 3J This type of elec- 
tricity transport does not require actual itinerant electron 
carrier movement, only local 7r-bond bonding energy ex- 
change. Third, analogous to the transverse wave in wa- 
ter, the locality of 7r-electrons can be viewed in parallel 
to the hydrogen-bond coupling among water molecules. 
Unless there exists a 7r-bond that is local and with weak 
bonding energy, such an energy relay mechanism would 
not be possible. The bonding/nonbonding mechanism of 
the 7r-bond loses no energy in the energy relay process 
when considering the initial and final states only, i.e., 
AH = thermodynamically, which is in agreement with 
the quantum tunneling description for the massless Dirac 
fermions. 7 



VII. CONCLUSIONS 

In summary, we have constructed a hybrid orbital 
model for Bi2Se3 based on crystal symmetry, atomic co- 
ordination, and the outer-shell electrons. The 7r-bond 
trimer for Se on each quintuple layer unit can be used to 
successfully explain the gapped semiconducting behavior 
for the bulk and the conducting surface. We believe this 
real space view from a chemical bond perspective should 



be complimentary to the understanding of topological in- 
sulator that has often been approached previously using 
the quantum spin Hall effect and band topology. The 
chemical bond approach could help us design additional 
topological insulators more intuitively and shift our focus 
to materials containing an open conjugated system that is 
constructed with 7r-bonds on the surface. It is clear that 
more experimental works are necessary to trace phenom- 
ena related to the existence and behaviors of 7r-bonds, 
especially more theoretical works to bridge the real and 
reciprocal views of topological insulators in parallel to 
graphene. 
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